%
% interp_matrix_v_6pt.m
%
%   Form the interpolation matrix for the vertical velocity (v)
%   which is stored at the horizontal edges -- uses a six point
%   interpolation scheme which is the adjoint of restriction
%
%   input -- Nx number of cells in x-direction, fine grid
%            Ny number of cells in the y-direction, fine grid
%   output
%
function Iv = interp_matrix_v_6pt(Nx,Ny);

% form two one-dimensional interpolation matricies
%
Ipx =  cc_nm_interp1_const(Nx);
Ipy = ec_interp1(Ny);

% assemble the matrix for the two-dimensional problem
%   note the x-y ordering in kron correpsonds to the ordering in meshgrid
%
Iv = kron(Ipx,Ipy);

